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ABSTRACT 

Using the fact that the null geodesies in NUT space lie on spatial cones, we 
consider the gravomagnetic lens effect on light rays passing a NUT deflector. 
We show that this effect changes the observed shape, size and orientation of a 
source. Compared to the Schwarzschild lens, there is an extra shear (a differ- 
ential twist around the lens axis) due to the gravomagnetic field which shears 
the shape of the source. Gravomagnetic monopoles can thus be recognized by 
the spirality that they produce in the lensing pattern. 

All the results obtained in this case (magnification factor, orientation 
of images, multiplicity of images, etc. ) depend on Q, the strength of the 
gravomagnetic monopole represented by NUT metric. One recovers the results 
of the usual Schwarzschild lens effect by putting this factor equal to zero. 

Key words: relativity - gravitation - (cosmology:) gravitational lensing. 

1 INTRODUCTION 

The usual gravitational lens effect is based on the bending of light rays passing a point mass 
M (Schwarzschild lens) in Schwarzschild spacetime . The static nature of the Schwarzschild 
metric implies that it can only produce gravoelectric fields, but in the case of NUT metric 
( sometimes called the generalized Schwarzschild metric ), the existence of the cross term 
ll dcf)dt n , shows that this space has also a gravomagnetic field B g (Demiansky & Newman 
1966). 

A simple example illustrates what we will analyse relativistically in the following sections. 
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Figure 1. Light ray passing a gravomagnetic monopole at large impact parameter b. 

Here we consider the deflection of an almost straight ray by the overall field of a stationary 
gravomagnetic monopole of strength Q and mass M. There are two fields of the following 
form: 

GMr , „ Qr 

E 9 = — aIld B g = ~^ 

and the corresponding force per unit mass on the ray is 

GMr Qv x r 

9 9 ' 

where v = c is the speed of the ray. The geometry is shown in Fig. 1. The ray passes the lens 
parallel to x-axis with a large impact parameter b. The normal gravitational force applies 
an inward radial force to the ray, but there is another force due to the gravomagnetic field 
along the z-direction. 

The deflection due to the gravomagnetic force can be found by calculating the impulse 

applied to the ray, assuming that it is undeflected: 

+00 +00 

AP 2 = / F z dt = Qb I <^^ = 2 _Q_ (1) 

The gravitational force deflects the ray in the x-y plane, but the gravomagnetic force deflects 
it out of the plane in the azimuthal direction so that the combined force twists the ray as 
the ray passes the monopole. 

2 NUT SPACE AND NULL GEODESICS 

The metric of the NUT space is given by (Kramer et al. 1980, Newman, Tamburino & Unti 
1963, Misner 1963) 

ds 2 = f(r)(dt - 2lcos6d<pf - f(r)- l dr 2 - (r 2 + l 2 )(d6 2 + sin 2 6d<j) 2 ), 
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Figure 2. The relation between r , j and L 



where 



f(r) = 1 



2(mr + I 2 ) 



(2) 



(r 2 + I 2 ) 



and 21 = Q is the strength of the gravomagnetic monopole represented by the NUT metric 
(7 is usually called the NUT factor or magnetic mass). Following the Landau and Lifshitz 
(1975) formulation of stationary spacetimes and using the standard variational method (i.e. 
S J ds 2 = 0), the authors (1996) have shown that all the geodesies in NUT space, including 
the null ones, lie on cones defined through the following equation: 



where L , r and e are the angular momentum, unit radial vector and the energy constant 
of the orbiting particle and j is a constant vector along the axis of the cone. Dotting equation 
(3) with r we find j • r = eQ showing that r lies on a cone; similarly L • j = L 2 = j 2 — e 2 Q 2 = 
const., so that L also moves around a cone. Later we will be concerned with the semi-angle 
of the cone or its complement X- From Fig. 2 one can see that 



L + eQr = j = const.S 



(3) 




where b = - is the impact parameter. In the limit Q -C b, i.e. small x we have 



Q 



t We use the natural units in which c = G = 1. 



© 0000 RAS, MNRAS 000, [j]-?? 



3 BENDING OF LIGHT RAYS IN NUT SPACE 

The equation of the null geodesies in NUT space is (Lynden-Bell & Nouri-Zonoz 1996) 

/(r)-V-f 2 )-^=0 , r = £, (4) 
where r is an affine parameter. 

Using the angle defined around the surface of the cone we have 

r 2 = L. (5) 
Now using (5) we can integrate (4) and obtain 



A0 



Ldr 



We perform the above integral in the limit where r 2 ^> I 2 + m 2 , i.e. far away from the event 
horizon (which is the case for the observed gravitational lenses) and where we have 

Ldr 

A0 



Changing to the variable u = - we have 

du 



r 

Umax 



n (e 2 — L 2 u 2 + 2mL 2 u i ) 2 



Calculating this integral to the first order in mu one can see that the deflection angle on the 
cone is 

4me 4m , . 

a = 2A0-7T = — — = — . (6) 

Lj 

which is the same as in the Schwarzschild case but with the difference that now b = - is the 

e 

impact parameter defined on the cone. 

4 GEOMETRY OF LENSING IN NUT SPACE 

The gravomagnetic lensing configuration is shown in Figs 3-5. A NUT lens is located at 
distances D d and D ds from an observer O and a source S respectively. The plane of the sky 
(or the lens plane) in which all the observations are made is the plane normal to Dd at N. 
For small deflection angle a, b <ti and D^, so from figure 3 we have 

P = ^r , & = t^- (7) 

© 0000 RAS, MNRAS 000, [|-?? 




A typical rim of the cone 
on which O and S lie 



Figure 3. A light ray (curved line) from S propagating in NUT space on a cone centered at N and passes the lens at distance 
b' , where b' 3> event horizon (To the first order in Xi b' — b, where b is the impact parameter). Note that S and O lie on the 
cone but not I, because OI is the straight light ray that the observer traces back to the image . 



Using these two angles one can see that 



Am b 

( p = 7T + a -/3-(3 s = 7i + — -— (1 + -=^-), (8a) 

and 

_j_ + + ft)] 

cosx cosx 

where <fi is the angle between generators and D^s around the surface of the cone and r\ 
is the projection of <fi onto the plane normal to the axis of the cone at its apex. 



5 LENS EQUATION 

The lens equation is basically the mapping of the source position to the image position on 
the lens plane. To find this mapping in the NUT case , we first need to find the angle between 
position vectors of the source and image (r and r') on the lens plane (see Fig. 5). 

As can easily be seen from the geometry in Fig. 5, the positions of the source and image 
on the lens plane are the intersections of this plane with lines OS and OI respectively (Points 
F and G). One should note that the point G is not on the cone but on the plane normal to 
the axis of the cone at its apex. 

From Fig. 5 one can see that the angle between the position vectors of the source and image 
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is the same as the angle between vectors r' and NM. Where NM is the component of D^s 
on the lens plane. Therefore, representing the vector NM by R, we have 



r'.R = r'.[D ds - (D ds .B d )t) d } = r'.D ds . (9a) 



On the other hand 



jD d cosx 



R = D ds sm6. (9c) 

where 9 is the angle between generators D d and D ds in the triangle ONS. Substituting 
equations (9b, c) in (9a) and using a spherical coordinate system centered at N one can show 
that 



(D d x D ds ).j 

D ds sxnVcosip = — — . 

jD d cosx 

But 



(D d x D ds ).j = \j\D d D ds cos 2 xsmr). 

Therefore 

cos 2 xsinr] 

cosip = — . (10) 

sinfc' 

Using the same spherical coordinates and equations (8a,b) one can show that 



CQ 3g = COsW" + a) ~ g|1 + (Dj/Dj -Vs.A. (11) 

cosx 

Equations (10) and (11) together determine ip the angle between the position vectors of 
the source and image in the lens plane. In what follows we will be concerned with the small 
angle (i.e. a,/3,x <C 1) approximation to the expressions (10) and (11) which are 

cos6 S -1 + 2 X 2 + \[a - (3(1 + ^)], (12a) 

and 
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Figure 4. Open, flattened cone and the light ray (dashed line) which is deflected at P passing the NUT lens, v = 2n(l — 
— j-) and a are the deficit and bending angles respectively. 

( L 2 +£ 2Q2 )2 - 

Lens plane 




1° 

Figure 5. Lens plane and the positions of source S, image I and the observer O. 



0(1 + Dd.) _ a - exL 

cos^ ^ — r - (126) 

[4*» + (a-/?(l + ^))V 
The following three limiting cases are worth considering: 

a) x = f° r which we get ip = for a > (3 and = 7r for a < (3, 

b) a = (3 = and }( < 1 for which we find that ijj — |, 
and finally, 

c) (3 ^> a & x from which we have ip = 2/j3(q/b) where (3 = (3(1 + j^)- 
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The above three results are all in agreement with what one expects from simple geometrical 
considerations. 

Having found the angle ip, the lens equation which maps the position of the source to that 
of the image, can be written as follows: 

A I r 

r = nAr k = — 

or 

{x = K,(x'cosib — y'sim/O 
(13) 
y = k(x simp + y'cos-0). 

where k is given in equation (17) and r and r' denote the position vectors of the source and 
image in the lens plane with the coordinates (x,y)(see Fig. 5). 

6 THE MAGNIFICATION FACTOR 

To see what happens to a bundle of rays in passing a NUT lens, we consider the magnifica- 
tion of an infinitesimal circular source centered at r; i.e. (Fig. 6) 

(X - xf + (Y - yf = 5x 2 + by 2 = s 2 s < 1. (14) 



Now by varying equations (13) and using (14) we have 

Sx 2 + 8y 2 = s 2 = (5k) 2 (x' 2 + y' 2 ) + k 2 (5x' 2 + 5y' 2 ) + /c 2 (<fy)V 2 + 

+2k(Sk)(x'Sx' + y'Sy') + 2k 2 (x 'by — y'5x')5i/j (15) 
On the other hand, using the geometry of lensing as shown in Fig. 5 for small angles, we have 

r' = (3D d S b, (16a) 

and 

r = lD d = D d X n A, ■ (166) 

where in finding the angle 7 we used equation (12a) and the fact that (fig. 5) 

1 D d = £D ds , 7 + f = 7T - 9. 
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Figure 6. Small circular source at r lensed to an ellipse at r' 



So 



and 



where 



r _ 7 _ [4 X 2 + (q-^] : 
r' /3 /? 



5 k 



,56 (2q/3 - 8k 2 - 2q 2 ) 
[4x 2 + (a-^)f ^ 



/8 = 0(l + -^4, 

Mis 



56 = Sr 



, (x'5x f + |/<V) 



and we have used the fact that for small angles a,f3 and x 

SP = P(j) , Sa = -a(j) , 6 X =-X(j)- 
Using equations (19) and (12) one can show after some manipulation that 

* = ^ -, A 

[4 X 2 + (a - ,3) 2 ] 6 

Substituting (17), (18), (19) and (20) in (15) and after a simple but long calculat 
finds that to the lowest order in angles and increments 

2 _ Aa{x'5x' + y'5y') 2 [A X 2 + {a - ffi 



b 2 (3 



ff 1 



\{8x'f + {8y'f]+ 



+ -^={x'8y' ~~ y'Sx')(x'8x' + y'dy'), 
Putting y' — 0, without loss of generality, the equation of the image becomes 



(A + B)X' 2 + BY' 2 + 2CX'Y' = 1, 
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where 

A=— C=i* 
/3s 2 ' (3s 2 ' 

(23) 

[4 X 2 + (g-^) 2 ] 

= =2 ) 

S 2 /3 

and we put 5x' = X' and 5y' = Y'. 

One can see that (22) is the equation of an ellipse in X'Y' coordinates. By definition the 

magnification '/x' is the ratio of the area of the image to the area of the unlensed source, so 
the area of the image 1 1 

^ ~ the area of the sourse ~ (A+A_)i ~ [(1 - ^)(1 - p - ' 
where As are the eigenvalues of the rotation matrix that transforms the ellipse into its 
principal axes (x'V).(Fig. 6) 

As it can easily be seen for \ = one recovers the known result of the Schwarzschild lens. 
The orientation of the ellipse is determined by the angle ( shown in figure 6 

2X 



tan( 



(a 2 + Ax 2 ) 12 - a 
or in terms of the strength of the gravomagnetic monopole Q 

tan( = ^—i • 

(4m 2 + Q 2 )2 -2m 

For Q = 0, as expected, we recover the result of the Schwarzschild lens which is C — f • So 
the presence of the NUT lens, not only changes the shape and size of the source but also 
its orientation and it is this last change which produces the spirality effect in the lensing 
pattern. For example for Qi = ^ and Q2 = m we get (j = 82.98 and Q2 = 76.71. Therefore 
the effect of the gravomagnetic field even in these extreme cases is quite small. This twist 
in the observed image would not be an unexpected effect if one looks at the null tetrad 
formulation of null geodesies in NUT space (see the appendix ). 



7 MULTIPLICITY OF IMAGES 

One can easily see that the limit \ = of the lens equation recovers the results of the 

Schwarzschild lens concerning the number of the images. In this section we will consider the 

number of the images for a given position of the source when x 0- 

Writing equation (17) and using the fact that to the lowest order in x, r' = b, we have 

r [g£ + (f ~ hAf\\ 
b bA 
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where A = D ^ + ^ d ■ Using the above equation one arrives at the following equation for b 

A 2 b' 2 + Bb' + AC = , b' = b 2 , (25) 



where 



B = -(8mA + r 2 A 2 ), 
C = 4Q 2 + 16m 2 . 



Solutions of (25) are 



where 



,4m r 2 , 1 r—.i , . 



A = r 4 A 4 + 16mA 3 r 2 - 16A 2 Q 2 . (27) 



We see that, in general, for a given r there can be two solutions for b (when A > or 

equivalently r 2 > r 2 ; see below). When A = 0, there is one solution for b which happens at 

9 Q 8m . Q 2 . i 8m 

In this case using (27) and (28) one can see that for Q = — > r = and therefore 

7 Am i 

» = (^-)-- 

This is the solution which corresponds to the Einstein rings in the Schwarzschild lensing. 
In figures 7 and 8 a horizontal line of point sources and its images for the two cases a < (3 
and a > (3 are shown and Figs 9-13 show randomly distributed small circular sources ( for 
both m = Q = 1 and m — , Q — 1 cases with —>■ oo ) in the plane of the sky and 
their lensing pattern. One should note that when m = 0, according to equatin (6), a = to 
the first order in — but by performing integral (5) to the lowest order in | (for m = 0) we 
can show that a — ~t$ (Appendix B). To have a better idea of lensing, some of the point 
sources as well as the circular ones and their corresponding images are labled. For a < (3 all 
the circles are magnified but for a > (3 ( Figs 11 and 13 ) all the circles are reduced. The 
spirality of the lensing pattern can be seen more clearly near to the center of the deflection. 
What the observer will actually see is the superposition of the two cases a < (3 and a > (3 
( Figs 10 and 12 ). 

As is clear from equation (26), for A < (or equivalently r 2 < r 2 ) there are no solutions 
for b, in other words an observer at a given position will not be able to see the source, if it 
lies inside a circle of radius less than tq. (point source a of figure 7 and circular source b of 
Fig. 9 are absent in the lensing pattern of the Figs 8 and 12 because of this and the same 
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reason applies to the incomplete images in Figs 10 - 13 ). This does not mean that there is a 
gap or a hole (in the case of a uniform distribution of point sources) in the observed image, 
because as is shown in Figs 10 and 12, the reduced images ( along with the images of the 
sources nearer to the lens which have smaller r by equation (28) ) will fill in the expected 
hole area. 

7.1 Part of the sky which can not be seen by a specified observer 

As we mentioned in previous section equation (25) has no solutions for r < r . This means 
that for an observer with a given distance from the lens a part of sky is unseeable. 
Expanding equation (28) in powers of £2 and keeping the first order term we have 

ro 2 



Q 2 



Am 

Figures 14 and 15 show r as a function of D ds , the distance between the source and the 
lens. 

7.2 Critical orbits and the dark area of NUT black hole 

To find the area at the NUT lens where due to the trapped null geodesies the observer sees 
(looking directly at NUT hole) a small black sphere we need to find the critical orbits which 
are specified by 



TT2 = ( 29 ) 



dr d 2 r 
d0 = ° #2 

Using (5) one can see that the above two equations are equivalent to the following two 
equations 



1 ,dr . 
-i 

d 2 r 



G(r) = -r^i^-) 2 = = (r 2 + I 2 ) 2 - (r 2 -I 2 - 2mr)b 2 = (30) 
~~ z + i z aw 



, = G'(r) = r 4 -l 4 - l 2 b 2 - mrb 2 = (31) 
dcp z 

Now substituting for r 4 in (31) from (30) we have 

r 2 (2l 2 - b 2 ) + 3mb 2 r + 2l 2 (b 2 + I 2 ) = 

Solving this for r we find 



3mb 2 + yWfr 4 + 8l 2 (b 2 + l 2 )(b 2 - 2l 2 ) 
T = 2{b 2 - 2l 2 ) (32) 

Now to find the impact parameter which is basically the radius of the dark area we need r 
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in the two limiting cases in which a) I <C m and b) m < L 
For (<mwe find from (32) 



2 I 2 
9 m 

substituting this into equation (31) and solving for b 2 we have 



3m(l + - — ) (33a) 



6 2 = 27m 2 [l + i(4)] (336) 
3 mr 

It is clear that these equations give the right answer for the Schwarzschild case i.e. for I = 

r = 3m , b = 3V3m (34) 

but because we took the limit in which I <C m, the above results do not give the exact value 
for the pure NUT case (m=0), for this we need to repeat the above procedure for the m <C / 
limit where from (32) we have 

_V3 i [l + ^(A D + -|-H)l (35a) 

b = 2V2 I [1 + D(m/ 1)] (356) 

in which D = 5 ^/| 8 - These equations are in agreement with the results of the pure NUT 
case which are 

r = V3 I , b = 2V2 1 (36) 

An interpolation formula between equations (33a,b) and (35a, b) which are good when either 
m or I is small is given by 

2 _ m 2 (9m 2 + 13Z 2 ) + 3/ 3 [/ + 2m(^D + ^ - §)] 

r - = pT^ 2 (37a) 

2 _ 9m 2 (3m 2 + Al 2 ) + 8/ 3 (Z + Dm) 

Equations (37) show that a circular orbit of radius r N is an unstable null geodesic and 
null geodesies arriving from infinity with an impact parameter less than b N cross the event 
horizon and get trapped. Therefore 6 at is the radius of the small black sphere the observer 
sees when he looks directly at the NUT black hole. 
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8 TIME DELAY BETWEEN IMAGES 

As we saw in the last section gravomagnetic lensing by NUT space can produce multiple 
images of a source. Different images correspond to different impact parameters and path 
lengths and therefore to different travel times. 

As in the usual lensing two different effects are responsible for this delay, namely the 
difference in path lengths and the difference in the gravitational potential traveresed by 
different light rays. In what follows we calculate the time delay between the arrival time for 
an image in the presence of the NUT lens and in its absence, taking into account both of 
the above contributions. 



8.1 Geometric time delay 

Geometric time delay by definition is 

(At) geo . = Al/c. (38) 

where A/ is the extra path length of the deflected light ray relative to an unperturbed null 
geodesic in the absence of the lens. Using Figs 3.3-3.5 and the geometrical relations defined 
there we see that 

Al = (D 1 + D 2 )-\T> ds --D d \, 

where 



D 2 = C0S ^-/ s) D ds = (! + -)[!- -(- - ^}D ds , 
cos% v 8 n 2 V 2 h ' J ' 



and 



\T> ds - D d | = (D ds + D d ) - ^f* [2x 2 + h a _ 

is the norm of the vector joining the source and the observer. 
Therefore 

A/ = l -D ds {P s a - (3 S 2 ) + ±D d {/3a - (3 2 ) + [2 X 2 + \{a- /3) 2 )]. 

Z l Ud + U ds £ 

Substituting for (3 and f3 s from equation (7) we have 

(A*W = Al/c = k\ 1 ^j V c? + 2( 7 ^^)x 2 ]- (39) 

C l D ds + L) d L> d + JJ ds 
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8.2 Potential delay (Shapiro's delay) 

As we showed in section 1.2.5, potential delay in NUT space can be calculated by using the 
Fermat generalized principle for stationary spacetimes (Landau & Lifishitz 1975, Schneider, 
Ehlers & Falco 1992) which reads; 

to = s(kf r- -i ■ Ta ^m-i ) = « = 1,2,3. (40) 



,cJ [goo' 1 + (A,^)]" 1 
where A a = ZL ^ L and dl is the Euclidean distance element along the light trajectory.^] 



900 

Equation (49) shows that goo 1 + A a ^£- can be thought of as the refractive index, n, of the 
spacetime under consideration, (in this case the NUT metric). Using this fact, the potential 
delay will have the following form 

, . . r dl r dl 

n 

or for the NUT space 

1 ^Observer 9m Hn 1 /-Observer 

(At) pot . = - [(1 - — y 1 - Qcos^dl -- dl 

C J Source V dl C ./Source 

where r , 9 and rj are the usual spherical coordinates. Using the facts that m<r, cosO = sin% 
and ri = we have 

' cosx 

1 /-Observer 9"W — 

(At) pot . = - —dl + ^(n + a-P)b, (41) 

C J Source T C 

where we used equation (8a) and the fact that tan% = x — §■ h can easily be shown that 

the first term to the first order in — is the potential delay for Schwarzschild metric, i.e. 

ip^ J(= 2™ h 4«. (42) 

C J Source T CO 

Therefore the total time dealy in NUT lensing is 

(At) to , = (At) 9eo . + (At) pot = + 2{-^-)\Y + 

C ZJJ ds L> d + U ds 

2M 4D d D ds x 2 , M s u 

Again we note that (43) reduces to the Schwarzschild lensing time delay when x = 0. 
Note that if the calculation were to be done in R-W spacetimes the D's would have been 
distances determined by the apparent size i.e. D d is evaluated at the time when light passes 
the lens ( redshift Z d ) and D ds is evaluated when light is emitted, (redshift Z ds ) 



Note that the calculations are done in the limit where r 2 2> Q 2 ,m 2 i.e far away from the horizon. 
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9 DISCUSSION 

We have studied the lensing of light rays passing a NUT deflector in the context of gravo- 
electromagnetism. We have shown that there is an extra shear due to the presence of the 
gravomagnetic field of the NUT space which shears the shape of the source. It was shown 
that this effect produces a spirality in the lensing pattern of a random distribution of circular 
sources located at infinity. This spirality effect about a NUT lens is very characteristic and 
is not displayed in normal lensing and the gravomagnetic lens due to a rotating object seen 
pole on does not show it because the twist of the ray as it approaches such a lens is cancelled 
by the opposite twist as it recedes. Thus the discovery of a spiral shear field about a lens 
would indicate, at least in principle, the presence of a grvomagnetic lens. 

The observability of gravomagnetic lenses is also discussed in the context of the con- 
straints on magnetic masses (Mueller & Perry 1986). It was shown by Misner (Misner 1963) 
that one can remove the string singularity in NUT space by using two different coordinate 
patches to cover the northern (0 < 9 < ~) and southern (| < 9 < tt) hemispheres at the 
cost of identifying the time coordinate with period 87ri]^ O ne can show that this periodicity 
is equivalent to the existence of closed timelike curves which in turn violates the causality. 
Mueller and Perry then discuss that since the Universe is not observed to be periodic on 
the scales less than 10 10 yr, any magnetic mass is presumably very large and at least 10 60 
Planck masses. This means that the characteristic scale size for such a mass would be the 
Hubble scale or larger. This large value for magnetic masses makes the lensing effect the 
most promising way to look for them. Although the expectation must be small, such effects 
should be looked for by those studying gravitational lenses. 
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APPENDIX A: NULL TETRADS AND THE GEOMETRY OF NULL 
GEODESICS 

Consider the family of null geodesies of a metric g ab . Let l a = ^f- denote the null tangent 
vectors to these geodesies, where u is an affine parameter along the geodesic. At each point 
one can define a complex, quasi-orthonormal null tetrad by the equation (Sachs 1961) 

g ab = 2&n b) -2rh (a m b \ 

where n a is another real null vector and rh a is the complex conjugate of m a . The above 
equations holds if and only if the quasi-orthonormality relations 

l a n a = -rh a m a = 1, 

l a l a = n a n a = m a m a = l a m a = n a m a = 0, 

hold. The vectors (l a , n a , m a , rh a ) define what is called a null tetrad. 

Using this null tetrad one can find some of the geometrical properties of the null geodesies, 
e.g. define the following two complex quantities (Hawking 1973) 

p = l a - b m a rh b , 
a = l a -b r m a m b . 

The real part of p measures the average rate of convergence of nearby null geodesies in space 

g ab and therefore is responsible for the magnification or reduction of the observed image. 

The imaginary part of p measures the twist or rate of rotation of neighbouring null geodesies 

i.e. this is the part responsible for the spirality in the lensing pattern shown in the text. One 

can easily calculate p for the Schwarzschild and NUT metrics: 

— r + il 

pNUT = ^Tp 
i 

PSchw. = 

r 

We note that the same correspondence between the two metrics also applies here i.e. 

Pschw. = hmpvf/T- 

The above relations show that in NUT case not only the size but also the orientation of the 
observed image changes ( Imp = r i l +l $ ) but in the Schwarzschild case there is no twist in the 
observed image because Imp = 0. The quantity a measures the rate of distortion or shear 
of the null geodesies. The effect of shear makes an originally small circular cross section of 
a bundle of rays become elliptical. 
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APPENDIX B: BENDING ANGLE FOR PURE NUT (M = 0) 

In this appendix we calculate the bending angle for a pure NUT case i.e. for m=0. Writing 
integral (5) for pure NUT we have 

Ldr 



r 2 (e 2 -£{l- 2 ^ 



where (1 — ^L) = e 2A | m=0 . Now expanding the integrand in powers of ^ and keeping the 



lowest order term we have 

A^ = b [ *a r 

J ( r 4_ & 2 r 2 + 2/2^2)5 

where b = L/e is the impact parameter. One can write the above integral in the following 
form 

rr m in=b (jf 

A<f> = b r (Bl) 

Joe ([ r 2_ a 2][ r 2_ d 2])2 

where 

a 2 = 2l2b ] b 2 (B2) 

b 2 /2-b/2Vb 2 -8P V ; 



d 2 = b 2 /2 - b/2Vb 2 - 8l 2 S 2l 2 (B3) 

Now due to the fact that the condition r m « n > a > d > is satisfied the above integral is an 
elliptic one of the first kind i.e. (equation 12 page 246 of Gradshteyn and Ryzhik 1980) 

= (b/a)F(u,t) (BA) 

where v = arcsin(a/r miri ) and t = d/a. Substituting (B2) and (B3) in (B4) we have 

A(j) = b/aF(ir/2, V2 l -) = b/a J -(T /2 [l - 2(/ 2 /6 2 )sin 2 #]- 1/2 ^ = tt/2 + vr/4(/ 2 /6 2 ) 

where we have used the fact that I ^ b. So the bending angle a is 

a = 2A0-7T = tt/2(/ 2 /6 2 ). 
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Figure 7. A horizontal line of point sources in the plane of the sky . 
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Figure 8. Corresponding lcnsed point sources for m = Q = 1 and for both a < /3 (primed ones) and a > /3 (double primed 
ones) cases . Note that the point source a can not be seen by the observer. Using a much denser line of point sources one will 
find that the two images form a continuous curve, i.e. there is a point between sources a and b whose two images coincide . 
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Figure 9. Randomly distributed circular sources in the plane of the sky . 
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Figure 10. Lensed circular sources in the plane of the sky for m=Q . The central part which is the image for the case a > j3 
is enlarged in the next figure . 
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Figure 12. Lensed circular sources in the plane of the sky for m=0 ( i.e pure NUT ). Note that the source b is completely 
absent in this case . The central part which is the image for the case a > /3 is enlarged in the next figure . 
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Figure 15. Three-dimensional version of fig. 14 for m = 1, Q = m/5, = 1. Note the difference of the scales on the Lens 
plane and along D^ s . 
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